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Abstract. We study completely asymmetric 2-cliannel exclusion processes in 1 
dimension. It describes a two-way traffic flow with cars moving in opposite directions. 
The interchannel interaction makes cars slow down in the vicinity of approaching cars 
in other lane. Particularly, we consider in detail the system with a finite density of 
cars on one lane and a single car on the other one. When the interchannel interaction 
reaches a critical value, traffic jam occurs, which turns out to be of first order phase 
transition. We derive exact expressions for the average velocities, the current, the 
density profile and the k- point density correlation functions. We also obtain the exact 
probability of two cars in one lane being distance R apart, provided there is a finite 
density of cars on the other lane, and show the two cars form a weakly bound state in 
the jammed phase. 



1. Introduction 

Low - dimensional systems out of equilibrium attract much attention recently 
An important class of such models is the one-dimensional (1-D) exclusion processes 
describing particles hopping independently with hard-core repulsion along 1-D lattice. 
Such systems provide good description of growth processes, traffic flow and queueing 
problems 0, etc. (see for the references up to 1995). The completely asymmetric 
exclusion process (ASEP) describing particles hopping only to the right with equal rate 
1 and hard core repulsion is perhaps the simplest and the best studied one In 
particular, for the periodic boundary condition, all configurations are equally likely in 
the steady state, and the average particle velocity in the infinite system is 



n being the density of particles. Janowsky and Lebowitz ||^ showed that a fixed blockage 
in the system that reduces the rate of hopping across it from 1 to r < 1 can produce 
global effects. Namely, for each fixed density n, there is a range of < r < ro where 
the system segregates in high- and low- density regions with sharp boundary, called 
shock, between them. Although some exact results were obtained 0, many quantities 
of interest e.g. steady state density profile, correlation functions, etc. were computed 
only numerically. 

In the present paper we derive all these quantities exactly, in closed form, for a 
slightly different model, guided by modelling the two-way traffic fiow problem. Namely, 
there are two 1-D chains on a ring, N sites each. One chain, or lane is occupied by 
cars and another with trucks (we call them differently just for notational simplicity), 
hopping in opposite directions with rate 1 and 7, respectively. Effective rate of hopping 
of the car (truck) reduces to l//3(7//3), when there is a truck (car) in front in another 
lane. Physically, 1/(3 is determined by narrowness of the road; it describes how much 




(1) 
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car/truck slows down seeing another truck/car approaching. ^ = case corresponds to 
the road completely blocked. 

For the case of a single truck in one lane and finite density of cars n in the other, we 
expect the similar type of behaviour with the blockage case . We expect the cars to 
pile up causing traffic jam, at a certain range of interlane interaction parameter (3. So it 
is, as Monte-Carlo simulations unambigiously show. Then, to study closely the nature of 
traffic jam phase transition, we impose the restriction that car and truck cannot occupy 
each neighbouring site i simultaneously. At this point, the model becomes exactly 
solvable by the matrix approach of Derrida et. al. P, ^. Using it, we compute the 
average velocities, the density profile, /c-point correlation functions exactly, for the finite 
chain and in the thermodynamic limit. Particularly, the traffic jam phase transition 
curve is given by simple formula (3cnt = ^/n/m the thermodynamic limit. 

Characterictics of the traffic jam phase transition are examined in detail. Finally, we 
consider the situation with two trucks/finite density of cars in the system. We observe 
that a weak bound state is formed between the two trucks. 

Comparison of our results with the original two-way model, the ASEP with single 
fixed blockage ||^, and exact Bethe ansatz solution by Schiitz |p for deterministic ASEP 
with blockage concludes the paper. 



2. Formulation of the 2-way traffic flow model 

We consider the following hard core exclusion process: there are two parallel 1-D chains 
on a ring, sites each, the first chain containing M cars and the second one containing 
K trucks. Cars (trucks) are hopping in opposite directions with rates 1 ( 7) respectively. 
The state of the system is characterized by the set of occupation numbers {tj}^^ of the 
first lane and {ai}f^i of the second lane. Xj = 1 if there is a car at site i and = if 
empty, and the same for the trucks, cTj = 1(0) if occupied (empty). The system evolves 
under the following stochastic dynamical rules: 

At each infinitesimal time interval dt, one pair of adjacent sites i,i + 1 at any of 
the two chains is selected at random for a possible exchange of states. The possible 
exchange processes together with their rates are listed below: 
cars are hopping to the right 

(Ti, Ti+i) = (1, 0) (0, 1) with rate | \ " ° , ^ . , , (2) 



if CTj+i = 1 (truck in front) 



trucks are hopping to the left 



(0,1)^(1,0) with rate <! { ^ " ' ° . • f ,^ (3) 

I ^ if r,- = 1 (car m front) 
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One sequence of 2N (total number of sites in two chains) selections constitutes one time 
step (or one Monte Carlo step, see Figures 1,2). 

The interlane interaction parameter (5 > 1 has a transparent physical meaning. It 
describes how much a vehicle slows down seeing another vehicle approaching, which in 
turn depends on narrowness of the road. (3 = 1 (no slowing down) corresponds to a 
highway with divider, and = corresponds to a narrow road completely blocked. 
If we let the system evolve for a long time, it reaches the steady state, independent of 
its initial configuration. We shall be interested in the steady state characteristics which 
depend only on macroscopic parameters (number of cars, number of trucks and total 
number of sites) and the rates (0,^. 

We studied the most practical characteristics of the model, the average velocities of 
cars (fear) and trucks (f truck) as a function of (1 — 1 //9) by Monte-Carlo (MC) simulations. 
The MC results for two different cases are shown in Figs. 1 and 2. Fig.l corresponds 
to the system with the density of cars and trucks being 0.3 and 0.4, respectively, and 
shows monotonic decrease of both velocities as j3 increases. Fig. 2 corresponds to the 
system with a single truck and many cars with density n = 0.3. In contrast to Fig. 1, 
we see that (fear) stays constant (equal to the average velocity (|I]) in the noniteracting 
system f = {1 — n), until the point 1 — J- ~ 0.8 is reached. For /3 > Pc, ("^car) rapidly 
drops. Simulations show that for (3 > (3c, the system segregates in two phases: the 
high density one in front of the truck (traffic jam) and the low density one behind the 
truck. Piling up of cars in front of the truck accounts for the decrease of average car 
velocity. Absence of sharp transition for a finite density of trucks (Fig.l) compared to 
Fig. 2 is due to the fact that a finite number of trucks in infinite systems produces the 
macroscopic jammed phase, while a finite density of trucks produces only microscopic 
jams which average out to give smooth behaviour. 

The segregated or traffic jam phase is well known as a shock phase or coexistence 
phase in 1-D asymmetric exclusion processes (ASEP). Schiitz |§] showed its existence in 
an exactly solvable deterministic model with a fixed blockage, and found various shock 
characteristics rigorously. As far as the probabilistic ASEP are concerned, Janowsky 
and Lebowitz ^ showed the existence of segregated phase in a probabilistic ASEP with 
a fixed blockage. The latter model is not solvable, and most results obtained in are 
therefore numerical. 

The characteristics of the shock are believed to be quite universal, qualitatively 
independent on details of stochastic process. That is why it is important to give exact 
solutions for some system with a shock. Here in this paper we propose a two-way traffic 
fiow problem (slightly modified, see below) with a single truck as an example of such a 
solvable system. Roughly, the single truck plays the role of blockage and ^ plays the 
role of transmission coefficient r in [^]. Following the traffic fiow formulation, we shall 
call the shock the 'traffic jam' and the segregated or coexistence phase the 'traffic jam 
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phase'. We shall find exactly the characteristics of the traffic jam and the traffic jam 
transition, including average velocities, density profiles, fc-point correlation functions, 
for finite chains and in the thermodynamic limit. To do so, one has to modify the 
original model to a solvable one. 



3. Modification of the original model to exactly solvable model 

We now modify the two-way traffic problem slightly. Here, we forbid the car and the 
truck to occupy the same site i in the neighbouring chains simultaneously. Then one 
can in fact describe the configuration by a single lane configuration {xj}^^; each site i 
is either occupied by a car Tj = 1 or truck Tj = 2 or empty Tj = 0. The allowed exchange 
processes is then modified from Eqs.(^ and to: 

(1,0) -> (0,1) with rate 1 

(0, 2) ^ (2, 0) with rate 7 (4) 
(1,2)^(2,1) with rate i 

Although the quantitative characteristics of the system do change after this modification, 
the qualitative characteristics do not (compare for instance the graphs for average car 
velocities in Figs. 1 and 3). 

The process (4) is the two-species ASEP solvable by the approach of Derrida et 
al [0. The process (4) and the one considered in 0] differ by replacement 2^0 
(interchange of trucks and empty spaces). 

The probability of a given steady state configuration is shown in Ref.[0 to be (up 
to normalization) the trace of a product 

Wconi{TiT2 . . . tn) = Tr (Xi X2 . . . X^) , (5) 

where 

{D if car at site i Tj = 1 
E if truck at site i = 2 (6) 
A if site i is empty Tj = 

are noncommuting matrices satisfying the following algebra: 
DE =D + E 

PDA = A (7) 

aAE = A; a = ^7 



Knowing the probabilistic measure, we can find the various averages of steady state. 
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4. Average velocities 

We shall consider the system having M number of cars and a single truck. 

Analogously to [0, define Y{N^M) as the probability of having the truck at site 
A, in a system with N sites and M cars. Define then Yd(A, M) as the probability 
of finding a car at site A — 1, provided the truck occupies the position A. Then, the 
average velocities of the cars and the truck are given by: 

, , _ 1 Yd{N, M) + {N-M- 1)Y{N - 1, M - 1) 
^^-■-Z p MY{N,M) 

_ 1 YpjN, M) + a (y (A, M) - YnjN, M)) 
l^truck;-^ F(A,M) ^> 

The quantities ^(A, M) and Ypi^N, M) are computed in Appendix and found to bef 
Y.(N, M) ^ -1, + "-^m M)) (10) 



1 



F(A, M) = Yn{N, M) + -^C^^_^ (11) 



a (3 



M 



where 



M 



I{N,M) = Y.P'C^:l, (12) 

k=l 

and Cl is the binomial coefficient. In Figs. 3 and 4 we plotted (fear) and (f truck), 
respectively, computed from exact formula, as a function of 1 — 1//9, for three densities 
n = 0.3, 0.5 and 0.7 at 7 = 1, in a system with A= 200 sites. Naturally the traffic jam 
transition point decreases as average density n increases. Behaviour of (fear) is similar 
to the one for the original two-way traffic model Fig.2. However now we can evaluate 
exact thermodynamic limit A, M — 00, n = M/N fixed, and find the exact transition 
point. We used the steepest descent method for computing the thermodynamic limits. 
The average velocities in the thermodynamic limit are given by (see Fig. 5): 

j 1-n iin(3<l 

(ftruck)= ? (l-n)(l-n,)+n(.+/3-„,) - ^ (^4) 



/3 



if n/3 > 1 



f More precisely, Y{N,AI) and Yd(^, Af) are probabilities up to normalization factor which is equal 
for all the terms entering Eqs.(^) and (||), see Appendix 
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Thus the transition point to the jammed state is given by a simple formula 

n/3crit = 1 (15) 

Note that the transition point does not depend on 7 — the free velocity of the truck. 
The average car velocity has a cusp at the transition point. The average velocity of 
cars before the transition nf3 < 1 is equal to the one in a system without truck (|l]). For 
n = 1, (vcar) = independently of because all sites are filled and cars cannot 
move. 1//? = is the case of complete blockage: both velocities are identically zero. 

To examine closely the nature of the traffic jam transition, we find the exact density 
profile and /c-point correlation functions in the next two sections. 



5. The density profile 

In this section, we obtain the exact density profile {n{x)) in a system with one truck 
and arbitrary M number of cars, in a chain of length N. We choose a reference frame 
in which truck is always at the position A^. It can be done because the weights of the 
steady state configurations depend only on the positions of the cars relative to the 
truck location, due to cyclic invariance of (|). The average density {n{x)) at distance x 
from the truck is equal to the probability of finding a car at site N — 1 — x; 

{n{x)) = ^conirN-i-.WM{nr,...r^ ^ 2) ^^^^ 

J2conf1^M{riT2. . .Tn-I 2) 

Sums run over all possible configurations having M cars, N — M — 1 empty spaces, with 
truck at the position A^, so 



wm = Tt\^{Y[X,) Ej (17) 

where Xi = D{A) if site i is occupied by a car (empty). The quantity (ITBI) is readily 
obtained using algebra (7) and found to be 



{n{x)) = {C^'zi - j^C^'-i + m - 1, M - 1) + 

(3''-\(3 - 1)I{N -x,M - x)e{M > x)]}/ {y{N, M)al3^') (18) 
where Y{N, M) and I{N, M) are given by Eqs. ( p!T]) and (1T2|), respectively, and 

e(i/>x) = | ^^f-^ (19) 

- ' I 0, otherwise ^ ' 

The expression (18) determines average density for any value of (5. Below we shall 
consider the cases before and after the traffic jam phase transition separately. 
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5.1. Low density phase; n[3 < 1. 

Before the transition the presence of truck affects the system only locally as is seen from 
Fig. 6, where the density profile is shown for = 400. The density, otherwise constant, 
locally increases only in a close vicinity of truck. 

Evaluating the formula (18) for large M, ^ 1; with ^ = n fixed, we find 

(nM) = n(l+ '°+/-y^-'" (n/3r) (20) 

Y 1 — 77, + an ) 

One sees that the local density disturbance decays exponentially at a finite length scale 



e=|ln(n/5)|-i (21) 

Therefore, the relative size of disturbed region vanishes as In principle, by common 
sense one would expect the existence of low density region right behind the truck. 
However it is absent in exact solution, as seen from Fig. 6. We do not have 
simple explanation for this fact. Analogous behaviour was observed in exactly solvable 
deterministic exclusion process with a fixed blockage [^]. 



5.2. Traffic jam phase; n(3 > 1 

In this region, cars pile up before the truck as seen from Fig. 7. Increase of interlane 
interaction /3 leads to increase of traffic jam length 

In thermodynamic limit M, A^ — > oo, ^ = the latter formula becomes exact, and the 
density profile becomes a step function: 

{■1 ^ < nl3-l 
1 \r '7' (23) 
■g otherwise. 

Note that the density in the low density region is equal to the critical density ricrit = 
independently of the average density n. The same behaviour is observed in 

In fact, in the jammed phase, the only way the truck can move is by the process 
(1,2) — ^ (2,1) (because the contribution of the processes (0,2) (2,0) becomes 
exponentially small in the large A^ limit). In the slow truck/ many cars problem 
treated in 0, the same is true when the rate for (2,0) — > (0,2) process as denoted 
by a in is zero. Thus our model in the jammed case is a special case of 0, up to 
exponentially small corrections. However only simple characteristics were investigated 
in M; correlation functions, as well as large A^ ^ 1 limits were not studied. 



More precisely, for large A^, M ^ 1, we find up to corrections of order N 2, 

= 1 - i(l - i) { 1 + erf (|-^) } + O (N-i) (24) 



y 



nf5-l ^ J2(3{l-n) 2 ,2 , 

with / = , A = -L and eii{y) = ^ e^* dt. 

P — I P — 1 VVTJO 



This shows that the shock interface extends over a region of width v N. More careful 
considerations however show that the real shock interface is sharp and extends over 
only 2 consecutive sites, and the apparent width of \/N is due to shock position 
fluctuations ( The shock position fluctuations of order were also observed in the 
probabilistic exclusion process with a fixed blockage 0). Indeed, the discrete version 
of the density gradient correlation, {An{xi) An{x2)) where An{x) = n{x + 1) — n{x) 
vanishes if \xi — X2I > 1. (It follows directly from Eqs. ( pT]) and (32)). This shows that 
the jammed phase is indeed segregated into two macroscopic regions — the low density 
one on the left with niow = ^ and the high density one on the right nhigh = 1- The fact 
that rihigh = 1 is due to the nature of the process we consider (see (4)): once cars pile 
up before the truck, the car-truck exchange processes do not create empty spaces. 

Finally, note that one can choose the difference between the average densities in 
the macroscopic regions in front and behind the truck dn = rihigh — ""-low as an order 
parameter, characterizing the traffic jam phase transition. With respect to this order 
parameter, the transition to the jammed phase is of the first order, as seen from Eg . (^Sf ) ; 



r ifn/5<i 

^"=1 l-i ifn(3>l ^25) 
5.3. The hydrodynamic approach 

The thermodynamic limit results Eq.(^) can as well be obtained from simple 
hydrodynamic arguments. Supposing that the segregated phase contains two 
macroscopic regions of length IN and (1 — /)A^, with average car velocities in these 
regions 1 — nhigh and 1 — niow, respectively (see Eq.(|T])), one can write down a set of 
equations. First, from the car conservation, 

raio„(l - /) + ^high^ = n 
and next, from the current conservation in the reference frame of the fixed truck, 

J = ^low(l - ^low + (^^truck)) = '^high(l - '^high + (?^truck)) (26) 

finally from the definition of the average car velocity 

n(fcar) = niow(l -/)(!- niow) + nhigh/(l - ^high) 
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Substituting the values (fear) and (f truck) from Eqs. (p!3D and (p^, and solving the above 
system of three equations, we obtain exactly the result Eq. (]23|) . 

The correctness of the hydrodynamic arguments in the thermodynamic limit is due 
to the fact that indeed cars behave like an ideal gas of interacting particles; correlations 
vanish in the thermodynamic limit as shown in the next section. 

One can easily compute the current flowing through the truck. The phase diagram 
"current versus ^ " (narrowness of the road) is given in Fig. 8. With fixed P, the current 
j increases with the density n, as 

j = n((t;truck) + (^^car)) = 

(1 a{l-n){l-np)+n{a + P-nP) \ 
""[P {l-n)il~nP)+n{a + P-nP) J ^ ^ 

until the critical density 

ncrit = ^ (28) 
is reached. After that, in the jammed phase, the current stays constant (see Eq.(PBD), 

imax = "-high(l - "-high + (^^truck)) = (29) 

for all densities ricrit ^ n < 1. 

Note that there is a single critical density value Eq. (pHj) , in constrast to the models 
with fixed blockage |^, where two critical densities exist, pcrit and pcrit = 1 — Pcrit- The 
reason is the following; the models considered in ^ have the particle-hole symmetry 
which is broken in the model we consider. 



6. The fc-point correlation functions 

Here we obtain the fc-point equal time correlation functions in the steady state, in exact 
and asymptotic forms, for a system with one truck and M cars. Analogously to Eq.(|TB|), 
one defines 

{n{xMx,) . . . n{x,)) = ^-^V" .^^-^ ^) (30) 

J2coniWM{m'2 ■ ■■T'N-l 2) 

with pj = N — 1 — Xj. 

Here we take Xi < X2 < ■ ■ ■ < x^j. Sums run over all possible configurations having M 
cars, — M — 1 empty spaces, with the truck at the position A^. 

f Equality of some argument values, say, xi = X2 simply lowers the order of the correlation function 
by 1 as seen from Eq.(|30|); {n{xi)n{xi)n{xz) . . . n{xk)) — {n{xi)n{x^) . . . n{xk)) 
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Calculation of Eq.(pO|) leads to the following surprising result: 

k 

{n{xi)n{x2) ■ ■ ■n{xk)) = {n{xk)) - fj{xk+i-j) (31) 

i=2 

The fc-point correlation function in fact splits into a sum of k terms, each one depending 
on a single argument! The exact form of fj{x) is given by 

X J{N -j,M- J, X - l)e(x > 2) + pC^'l^L,]}e{x > 1)]/ {y{N, M)af3'') (32) 
where Q{x>y) and Y{N^M) are given by Eg. (p!9D , (pT]), respectively, and 

min {M,x) 

JiN,M,x)= P'C^'IL^■ 

i=l 

We shall show that in the limit of large N,M ^ 1, fj{x) is given by a remarkably 
simple formula 

fj{x) = K^^^ (1 — {n{x))) ; with k = min ^n, — ^ (33) 

Indeed let us consider the connected two-point correlation function, 

{n{xi)n{x2))c = {n{xi)n{x2)) - {n{xi)){n{x2)) = 

(n(x2))(l-(n(xi)))-/2(xi) (34) 

according to Eg . (^T]) . As f2{xi) does not depend on X2, one can choose any convenient 
X2. Take the point X2 infinitely far apart from Xi, X2 ^ Xi, so that the correlations 
between them vanish {n{xi)n{x2))c = 0. Then 

a) before the transition nj3 < 1, we have for N, M ^ 1, using Eg . (pO]) and (34), 

/^(xi) = n{l - {n{xi))) 

b) after the transition n/S > 1 , using Eg. ( ^4]) and imposing in addition X2 ^ Nl we 
obtain {n{x2)) = ^, and 

/2(xi) = ^(l-(n(xi))) 

which proves the formula Eg.(^) for j = 2. Recursively, one obtains the asymptotic 
behaviour of other functions /^{x), . . . , fk{x) from (|3l[) . The formula (|33D can also be 
obtained directly from (32). 

Finally, the /c-point correlation function, connected part, is given from Egs.(pl|) ([33|) 

as 

k 

(ra(xi) . . . n{xk))c = {n{xk)) - (1 - {n{xk+i-j))) - 

i=2 
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. . . {n(xk)), where k = min ( ^, -5 1 (35) 




for N, M 1, xi < X2 < ■ ■ ■ < Xk and {n{x)) is given by ( |20D and (|2^ ) for < 1 and 
n(3 > 1 respectively. Thus the A;-point correlation function for A^, M ^ 1 is determined 
completely by the 1-point correlation functions . 

As an example, consider the two-point correlation function, 

{n{xi)n{x2))c = i{n'{x2)) - k)(1 - (^(a^i))) 

Before the transition, at the low density phase nP < 1, n = n, (n(x2)) is given by (j20|), 
and (n(x2)) — k ~ {n[3Y'^ . So the correlation function decays exponentially with a 
length scale 

e=|ln(n/5)|-i 

Thus, in the low density phase, the two-point correlation function is nonzero only in a 
close vicinity of truck xi < X2 ~ ^. For X2 ^ {n{xi)n{x2))c = 0. Thus, in the whole 
region ^ -C X2 < A^, cars do not feel any correlations between each other and behave 
like an ideal gas of particles. 

In the traffic jam phase, njS > 1, the two-point correlation function 

{n{xi)n{x2))c = i{nix2)) - -^)(1 - {n{xi))) (36) 

stays nonzero, only if both Xi and X2 {xi < X2) are in the region Xi,X2 € [Nl — 
A^/N,NI + Av^] as is seen from (p^. Otherwise, either {n(x2)) ~ ^, or {n(xi)) ~ 1 
and the correlation function {n{xi)n{x2))c vanishes. 

Again, one can say that the jammed phase indeed has a phase separation: a) solid- 
like phase with density 1. b) ideal gas phase (no correlations between the particles-cars) 
of density ^. This explains why the simple hydrodynamic approach (see section 
leads to the correct results in the thermodynamic limit. 



7. The bound state between two trucks 



Here we shall consider the system having two trucks, and M cars, and determine 
the probability fl{R) of two trucks being at distance R apart. This probability is 
proportional to 

Q{R) ~ ^m{tiT2 . . . tm-r-2 2 tn-r ■ ■ ■ tm-1 2), (37) 

conf 

with the sum running over all possible configurations having M cars and two trucks at 
positions (A^ — i? — 1) and A^. So = corresponds to two trucks being next to each 
other. Due to the periodic boundary condition, < < ^^y^- 
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The exact expression for Q{R) at finite N is unwieldy and we shall not present it 
here (Typical behaviour of Q{R) are shown in Fig.9 for = 200, n = 0.3, and /5 = 3, 5) 
Instead we shall write down its asymptotics in each phase. 

i) 71/3 < 1 : In the thermodynamic limit, Q{R) reduces to the following; 

^{R) - 1 + + (38) 

[1 — n + any 

It is maximal for = and decays exponentially with the same length scale ^ = 
I ln(ra/3)|~^, as before. As ^ does not depend on A^, the fraction of space with nonzero 
correlations between the trucks vanishes as 1/A^. Thus two trucks are asymptotically 
free. The asymptotic freedom of trucks accounts for the fact, that the phase transition to 
the jammed phase takes place at the same critical density, ricrit/? = 1- These arguments 
can be extended to any finite number of trucks in the infinite system N oo. 

ii) Jammed phase n(3 > 1: In the thermodynamic limit, we find the following result 
for VL{R) from the exact formula: 

Vl{R) linearly drops with the distance R in the region 

1 nS — 1 

0<R<Nro, ro = min(/, 1 - /) < -, / ^ 



2' (3-1' 

and then stays constant Q{R) = const, Nvq < R < ^y^, see Fig. 9(b). The relative 
ratios are 

njNr,) (a-l)(/3-l) ro , 

for ro = 1 - /, = 1 40 

iZ(0) ap 1 — tq 

One can interprete this as the two trucks forming a weak bound state in the traffic 
jam phase. The probability Q{R) was studied also in paper 0, for a system with two 
second class/many first class particles, both hopping in the same direction, where it 
shows a power law decay f2(-R) ~ in the uniform background of the first class 

particles. This is in marked contrast to our result Eq. (|38| ) showing exponential decay in 
the uniform low density phase. Note however that the asymptotic Eq.(p8D was obtained 
in supposition a ^ 1, /3 7^ 1, while the results of the paper are derived for a = 1, /5 = 1 
case. 

It is interesting to analyze the average distance between two trucks. Consider the 
case n < 0.5 first. Analysis of (^) and (^) shows that in the thermodynamic limit the 
relative distance between the trucks = 1/4 in the low density phase njS < 1 and 
then drops monotonically as function of /5 from = 1/4 to = n/3at/5 = cxo 
(complete blockage). At the same time the length / = {n(3 — — 1) increases from 

/ = at the transition point /?crit = l/n to l = n at (3 = 00. We recall the reader that 
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/ = '^^Zi is the total length of traffic jam in a system with 1 truck in the thermodynamic 
limit, see (|2^). One can argue in a different way that the total length of the traffic jam 
is independent on the number of trucks, as long as it remains finite, (e.g. by using the 
hydrodynamic approach, see subsection ^]3| ). That means in the early stage of traffic 
jam phase (small j3) > / and two separate traffic jams in front of the two trucks 
are formed, of lengths li and I2, h + h = I separated by the low density regions. As 
/3 increases, eventually < /, and the two separate jams merge into a single one of 
length /. 

8. Summary 

We have formulated the two-lane traffic model and showed that it has the transition 
from the low density phase to the segregated (traffic jam) phase. Modifying the model 
to an exactly solvable one, we studied in detail the characteristics of the solvable model, 
which we believe to be qualitatively correct for the original one. The solvable model in 
the traffic jam phase in ^ 1 limit is a special case of the 2-species model considered 
by Derrida see discussion after Eq.(^). However our angle of view is different and 
most results obtained in sections 4-7 are new. We have obtained exact expressions of 
the current j, the average density profile, and the fc-point correlation functions, for finite 
chain, and in the large N limit, for a single truck. We have also studied the two trucks 
case and observed that a weakly bound state is formed between them in the traffic jam 
phase. Generally, the truck slowing down the car movement, can be thought of as a 
sort of moving blockage. Qualitatively our results for the density profile, the current j 
phase diagram, two-point correlation functions agree with those obtained in and in 
part with those in ||^, describing a fixed blockage. However the last two systems possess 
the particle-hole symmetry and therefore the uniform high density phase, related to the 
low density one by this symmetry. In our case, the particle-hole symmetry is broken for 
both our original two-lane model and the modified solvable one. This accounts for the 
absence of the uniform high density phase in our model, see phase diagram Fig. 8. 
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Appendix A. Computation of Y{N, M) and Yd{N, M) 

The probability of finding the truck at the site is given up to normalization by 

y (iV, M) = Y, WM{nT2 . . . TM-l 2) (Al) 
conf 

where this sum is over all possible configurations having M cars, N — M — 1 empty 
spaces, with the truck at the position A^. Split the above sum in two terms as 

12) + J2wUriT2...TN-2 2) (A2) 

conf conf 

The first ( second ) term corresponds to a car (empty space) being at site (A^ — 1). The 
second term can be written as 

wUnT2 . . . TN-2 2)= Ti{CAE) = -TiiCA) (A3) 

a 

using the algebra (7). Here, C corresponds to arbitrary configuration of length N — 2 
having M cars and Q = N — M — 2 empty spaces. Generally, 

with mi + 1712 + . . . + rrik = M, qi + q2 + . . . + Qk = Q, Qk > ^■ 
According to (7), 

T^^CA) = Tr ( ^^,^^1^...^^, ^"+^-+-+^-) = ^Tr (A4) 

and 

wUnT2 . . . TN-2 2) = ^Tr (A5) 



The first term in (|A2|) for some specific configuration is 

WM(rir2 . . . tn-2 1 2) = Tr (D'''^ A'^W'' A"' . . . D'''^-' A'^'^-W"^ E) (A6) 

where mi + m2 + . . . + = M, gi + g2 + • • • + Qk-i = Q, ^k > 1- We have the following 
recursive relation: 

= AD'^E = AD'^-^DE) = AD"'-\D + E) = AD"" + (A7) 

Using the last expression recursively, one finds 

f^ = AY,D' + /o; fo = AE = -A; (A8) 

Substituting the value of fm into (|A^), using (7), we obtain 

< = ^^-^Tr (A9) 
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where 

7„. = a E + /3- = -__ + " + (AlO) 

The same common factor Tr cancels from all formula for averages starting 

from (H,^ and on, as it enters to both the denominator and numerator. Below we shall 
set Ti[A9^ = 1 for simplicity. 

Using ( |A5| ) and and some combinatorics to count the number of configurations, 
the sum (X2) then reads 

1 ^ 1 

M 

Finally, substituting ( |A10| ) and performing the summation CfJZm^2 = Cffzi, one 

m=l 

obtains Eqs.([Tl|) and 
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Figure captions 

Fig. 1. MC simulation result of the two- lane traffic flow model. Shown are the "average 
velocities versus interlane interaction r = 1 — ^" for = 200, ricars = 0.3, ntmcks = 0.4 
and 7 = 1 f. Initial conflguration of the system is random. We equilibriate the system 
for 2000 MC step intervals, and collect data at 2000 MC step intervals and average over 
100 different histories. 

Fig. 2. The same as Fig. 1 but for Ncars = M = 60 and Ntmcks = 1- The point Vc ~ 0.8 
is the approximate traffic jam transition point. 

Fig. 3. Average velocities of the cars computed from Eq. (P) in a system of 200 sites, 
for different densities n = 0.3, 0.5, 0.7. 

Fig. 4. The same as Fig. 3 for the the truck velocities Eq. (Qj. 

Fig. 5. Exact car and truck velocities in the thermodynamic limit Eqs.(|l^) and (Q), 
for n = 0.3. (ftruck) is given for 7 = 1 and 7 = 1.5 (thin line). 

Fig. 6. Density proflle before traffic jam, in a chain of 400 sites computed from formula 
(18). n = ^ = 0.3, P = 3 (/?crit = 3.3333 . . .). The single truck is located at right end. 

Fig. 7. Density proflle in the traffic jam phase, in a chain of 400 sites computed from 
formula (18) for n = ^ = 0.3,. and for P = 5(a) and /? = 15(6). {Pent = 3.3333. . .). 
The total length of traffic jam is Ljam ~ ^^^Er-^- Asymptotic values of densities in a 
high (low) density regions are ^high = 1, ""-low = ^- Truck is located at the right end. 

Fig. 8. Phase diagram of the solvable model in "current versus ^ " plane. The curve 
j = ^ corresponds to jammed phase niow = ^; ^high = 1- The area below corresponds 
to uniform low density phase niow = "'^high < ""-crit- There is no uniform high density 
phase like in |^, H because the particle-hole symmetry is broken. 

Fig. 9. Probability of flnding two trucks at a distance R apart, before (a) and after (b) 
the phase transition, computed from the exact formula, for = 200, n = 0.3. The 
thin line in (b) corresponds to the thermodynamic limit. The value of /3 is 3 and 5 for 
(a) and (b) respectively. 
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